Aliquots
These are numerous so probably require a database or other tiling strategy to serve (vs native R).
Renewable Energy
Renewable Energy - GIS Data | BOEM
Renewable Energy Leases and Planning Areas
Boundaries of renewable energy lease areas, wind planning areas, and marine hydrokinetic planning areas.
Code
# libraries ----
librarian:: shelf (
dplyr, glue, janitor, here, knitr, mapview, rmapshaper, sf,
quiet = T)
Code
# paths ----
dir_data <- "/Users/bbest/My Drive/projects/msens/data"
wind_gdb <- glue ("{dir_data}/raw/boem.gov/BOEMWindLayers_4Download.gdb" )
(d_wind_lyrs <- st_layers (wind_gdb, do_count = T))
Driver: OpenFileGDB
Available layers:
layer_name geometry_type features fields
1 BOEM_MHKLeasesandPlanningAreas Multi Polygon 17 21
2 BOEM_Wind_Planning_Area_Outlines_8_30_2023 Multi Polygon 23 8
3 BOEM_Wind_Leases_9_27_2023 Multi Polygon 3880 18
4 Wind_Lease_Outlines_9_27_2023 Multi Polygon 41 15
5 BOEM_Wind_Planning_Areas_10_06_2023 Multi Polygon 8685 11
crs_name
1 WGS 84
2 WGS 84
3 WGS 84
4 WGS 84
5 WGS 84
Code
wind_lyrs <- d_wind_lyrs$ name
ply_wind_lse <- read_sf (wind_gdb, wind_lyrs[3 ]) # BOEM_Wind_Leases_9_27_2023
ply_wind_pln <- read_sf (wind_gdb, wind_lyrs[5 ]) # BOEM_Wind_Planning_Areas_10_06_2023
ply_wind_lse_s05 <- ply_wind_lse |>
ms_simplify (keep= 0.05 , keep_shapes= T)
ply_wind_pln_s05 <- ply_wind_pln |>
ms_simplify (keep= 0.05 , keep_shapes= T)
mapView (
ply_wind_lse_s05,
layer.name = "Wind Leases" ,
col.regions = "red" ) +
mapview (
ply_wind_pln_s05,
layer.name = "Wind Planning Areas" ,
col.regions = "green" )